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Abstract. It is argued that an important fraction of PNe present large 
temperature variations that are not due to observational errors nor to 
incomplete atomic physics. Seven possible causes for these variations are 
reviewed, one of them is presented for the first time. 



1. Overview 

Photoionization models for chemically homogeneous gaseous nebulae of constant 
density predict an almost constant temperature and consequently very often 
observers assume a constant temperature to determine chemical abundances. 
There is growing observational evidence that indicates the presence of large 
temperature variations in gaseous nebulae, in contradiction with the models 
mentioned above. Consequently the relevance of temperature variations in the 
abundance determinations has to be considered in detail. 

In section 2 we review the observational evidence in favor of large temper- 
ature variations. In section 3 we review different possible explanations for the 
observed variations. In section 4 we discuss the implications of the temperature 
variations in the abundance determinations, these implications depend on the 
source of the temperature variations. 

Recent reviews on the temperature structure of gaseous nebulae are those 
of Peimbert, M. (1995, 2002), Esteban (2002), Stasihska (2002), and Liu (2002a, 
2002b). 



2. Temperature Variations 
2.1. Definitions 

The average temperature, Tq, and the mean square temperature fluctuation, , 
are given by 

!Uv)N,{v)N,{v)dV 
j N,{v)N,{v)dV ' 

and 

2 _ !{T,-T^fN,NdV 

respectively, where A^e and Ni are the electron and the ion densities of the 
observed emission line and V is the observed volume (Peimbert 1967). 
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To determine Tq and t"^ we need two different methods to derive Tg. one 
that weighs preferentially the high temperature regions and one that weighs 
preferentially the low temperature regions. For example the temperature derived 
from the ratio of the [O III] AA 4363, 5007 lines, 1(4363/5007) ) the temperature 
derived from the ratio of the Balmer continuum to I (HP), T'(Ba,c/if/3)5 that are 
given by 

r 1 /nnann \ 1 
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^(4363/5007) — ^( 

and 

T(Bac/H/3) = r(Bac) = ro(l - 1.70*^), (4) 

respectively. It is also possible to use the intensity ratio of a collisionally excited 
line of an element p + 1 times ionized to a recombination line of the same element 
p times ionized, this ratio is independent of the element abundance and in the 
low density limit depends only on the electron temperature. Two examples of 
this method are: 

^(O Ilrec/O lllcoll) = ^(4649/5007) = /l(^0)^^) (5) 

and 

^(C Ilrec/C lllcoll) = ^(4267/1909) = /2(7o,t ), (6) 

by combining these ratios with a temperature determined from the ratio of two 
collisionally excited lines, like ^(4363/5007), it is also possible to derive To and t^. 

From the ratio of two He I recombination lines it is possible to obtain 
relationships of the type: 

^(Hc ^ fijC^O, t^, ^6,^3889), (7) 

and from a set of 5 or more helium lines (which should include AA 3889, 6678, 
and 7065) it is possible to obtain Tq, , Ne, and r3889 (see Peimbert, Peimbert, 
& Ruiz 2000). Sometimes the errors in these determinations are quite large, in 
such case it is possible to combine the information from the helium lines with the 
one from T'(43g3/5oo7) to obtain better determinations of Tq and (see Peimbert, 
Peimbert, &; Luridiana 2002). 

2.2. Errors in the temperature determinations 

Different types of errors are present in the determinations of the electron tem- 
peratures, and in some cases these errors have led to erroneous temperature 
determinations. Before accepting that temperature variations are present in a 
given object, it is necessary to rule out any possible sources of error. 

Rola Sz Stasinska (1994) and Kola & Pclat (1994) have argued that there 
is a systematic overestimation of the line intensities of lines with low signal to 
noise ratios that leads to the overestimation of the ^(4267/1909) temperatures. 
This effect is present in many objects and has to be taken into consideration, 
but for many well observed objects it is not the main cause of the derived 
large T'(43g3/5oo7) - ^'(4267/1909) differences (e.g., Peimbert, Torres-Peimbert, & 
Luridiana 1995). 

The discrepancy between the C abundances derived from the C II A 4267 
recombination line and those derived from the collisionally excited 1906+1909 
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C III lines has been known for at least 20 years and the idea that the large 
intensity of the C II A 4267 line was due to an incomplete knowledge of the 
relevant atomic physics has been mentioned on several occasions. Nevertheless 
similar excesses in the intensities of recombination lines are observed in different 
multiplets of the same ion and in other elements (Liu et al. 1995; Mathis, 
Torres-Peimbert, & Peimbert 1998). 

Viegas Sz Clegg (1994) suggested that the differences between ^(4049/5007) 
and T(Bac) could be due to density fluctuations diminishing the 1(4959) and 
/(5007) values by collisional de-excitation yielding a spuriously high ^(4549/5007) 
value. This is indeed the case for the planetary nebula M2-29, where from 
ground based observations of the whole object a ^(4353/5007) of 24,000 K and 
a Nf, = 2000 cm~^ are derived, while from HST observations it is found that 
it has at least two density components (one with Ne ~ 10^ cm~^ and another 
with Ne ~ 10^ cm~^, the in the low density component is of 8800 K); the 
high temperature derived from the ground observations is spurious and yields 
0/H= 5.1 X 10~^, a factor of ten smaller than the value derived adopting a 
two density component model (Pefia, Torres-Peimbert, & Ruiz 1991; Torres- 
Peimbert et al. 2002). A similar result has been found for Mz 3 by Zhang &; 
Liu (2002). 

2.3. Observations of Temperature Variations 

Prom uniform density chemically homogeneous photoionization models it has 
been found that is in the 0.000 to 0.025 range, with typical values around 
0.005 (e.g., Gruenwald & Viegas 1995). 

From observations and by combining equation 3 with any of the equations 
4 to 7 it is possible to determine and Tq. These combinations indicate values 
of in the 0.00 to 0.15 range with typical values around 0.04. The errors in the 
best determinations given in the literature are typically in the 0.01 to 0.02 
range. 

Some of the best determinations of and Tq are: from equations 3 and 4 

those by Peimbert (1971), Liu & Danzigcr (1993), Liu ct al. (1995, 2000, 2001); 
from equations 3 and 5 those by Peimbert, Storey, & Torres-Peimbert (1993), 
Liu et al. (1995, 2000, 2001); from equations 3 and 6 those by Rola &; Stasinska 
(1994), Peimbert, et al. (1995); and from equations 3 and 7 those by Peimbert 
et al. (1995, 2000, 2002). 

Liu et al.(2001, see also Peimbert, A. 2002), presented a good correlation 
between Log{0~^^ /'5i~^)RL/{0^^ /}i^)cEL and T(43g3/5oo7)-T'(Bac) and mention 
that this correlation strongly supports the idea that the temperature variations 
are real. In Table 1 we present the ^(4303/5007) and r(Bac) temperatures de- 
rived by Liu et al. (2000, 2001) for three of the objects that show very large 
temperature variations. From these values we have computed T'(4g49/5oo7) and 
from equations 3 and 4 we have derived t^(3,4), and from equations 3 and 5 
we have derived t^(3,5). Prom this table it can be seen that t^(3,4) and t^(3,5) 
are not only proportional to each other but that they are in good agreement 
with each other. This implies that most of the temperature differences among 
r(Bac), T(4363/5oo7)) and T'(4649/5oo7) are due to real temperature variations in 
the nebulae, thus an explanation for the origin of these variations has to be 
sought. 



4 



Torres-Peimbert &: Peimbert 



Table 1 

Mean Square Temperature Fluctuation Values 



Ml-42 9220 3560 5550 0.123 0.120 

M2-36 8380 5900 6025 0.053 0.078 

NGC 6153 9140 6080 6220 0.065 0.097 



3. What Causes the Temperature Variations? 

In what follows we present seven different mechanisms that can produce large 
temperature variations. It is clear that all of them should be considered. Their 
relative importance probably varies for each object, and in many cases, the ones 
that show t"^ < 0.01 values, all of them are negligible. 

3.1. Shadowed Regions 

To explain the intensities of lines of low degree of ionization Mathis (1976) 
suggested the existence of regions shadowed from the central stars, these regions 
would be ionized by diffuse radiation and would have a considerably smaller 
temperature than those regions ionized directly from the central stars of H II 
regions. 

3.2. Chemical Inhomogeneities 

Some PNe, like A30 and A78, present important chemical inhomogeneities; 
therefore the idea that temperature variations can be due to various degrees 
of chemical inhomogeneities has been proposed. Chemically inhomogeneous 
photoionization models have been presented by different authors (e.g., Torres- 
Peimbert, Peimbert, &: Peiia 1990; Kingdon &; Ferland 1998; Pequignot et al. 
2002a, 2002b; Tylenda 2002). 

3.3. Density Variations 

Extreme density variations are present in most PNe, as can be seen from optical 
images; these clumps, shells, and filaments, can produce temperature variations 
by themselves under the assumption of photoionization heating only. Mihalszky 
&; Ferland (1983) have presented chemically homogeneous photoionization mod- 
els, including density variations in the 10^ to lO^-^cm"^ range, and find values 
of of about 0.005, which could explain only a small fraction of the largest 
discrepancies observed. 

3.4. Deposition of Mechanical Energy 

Shocks have been proposed as a mechanism to account for temperature fluctua- 
tions (Peimbert, Sarmiento, &: Fierro 1991; Peimbert et al. 1995 and references 
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therein). Peimbert et al. (1995) have found that those objects that show higher 
velocity dispersions are the ones that show higher temperature variations, sup- 
porting the presence of shocks. In a PN a shock front is expected to be formed 
at the interface between a fast stellar wind from the central star and the slowly 
expanding circumstellar envelope. The fast stellar wind would carve out a hot 
bubble in the interior of the nebula that would emit primarily in X-rays. Clear 
evidence for such hot bubbles has been found recently for a small number of 
PNc from X-ray imaging (Kastner et al. 2000; Kastner, Vrtilek, & Soker 2001; 
Chu ct al. 2001). 

Mellema (1997) finds strong temperature variations from hydrodynamical 
simulations of bipolar PNe. The dense gas near the equator has temperatures 
around 9000 K while the gas near the poles has temperatures between 20,000 
and 50,000 K. The reasons are two: at the poles the nebula has a higher shock 
velocity, leading to more heating; and a lower density, leading to less efficient 
cooling. Bohigas (2002) finds that in the outer regions of 9 PNe of Type I 
the r(0++)/r(N+) is large and the Ra/[S II] line ratio is small, a typical 
combination in shock excited plasmas. 

3.5. Deposition of Magnetic Energy 

Ferland, in the meeting on Ionized Gaseous Nebulae (Mexico City 2000), sug- 
gested that magnetic reconnection can provide localized temperature variations 
in ionized plasmas. According to Garcia-Segura, Lopez, &; Franco (2001) the 
multiple, regularly spaced concentric shells around some PNe could be due to 
the effects a of a solar-like magnetic cycle, with periodic polarity inversion, in 
the slow wind of an asymptotic giant branch (AGB) star. Presumably these 
shells of alternating polarity could give rise to magnetic reconnection processes 
once that they arc compressed in the formed PN (i.e. in the swcpt-up shell). 
Magnetic fields of the order of milligauss have been measured in the torus sur- 
rounding K3-35 (Miranda et al. 2002) ; also the presence of a magnetic field has 
been inferred in the planetary nebula 0110.9+1.3 (Zijlstra et al. 1989) and in 
the pre-planetary nebula IRAS17150-3224 (Hu et al. 1993). 

3.6. Dust Heating 

Stasiriska & Szczerba (2001) have analyzed the effects of photoelectric heating 
by dust grains in photoionization models of PNe. They have shown that this 
process is particularly important if filamentary PNe contain a population of 
small grains. The temperature structure of such dusty and filamentary nebulae 
would produce large temperature variations. 

3.7. Decrease of the Ionizing Flux 

The ionizing flux from the central stars of PNe is not constant. As a first 
approximation the ionization increases, then reaches a maximum and finally 
decreases. During the decrease and in those directions where the nebula is 
density bounded there would be ionized regions that would become isolated from 
the stellar ionizing photons, these regions would cool first and would recombine 
later (e.g., Spitzer 1978). A typical recombination time is ^ lO^/iVg years while 
the cooling time is an order of magnitude smaller, therefore cool regions with a 
relatively high degree of ionization are expected to be present. 
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This effect might be important for PNe produced by massive stars, because 
the central stars evolve relatively fast. One possible example of this mechanism 
might be provided by N66 in the LMC. In this object the ionizing flux has 
diminished by a factor of two in the last ten years (Pena et al. 1997; Pena, 
Hamann, & Ruiz 2002) and the RL/CEL abundance ratio is of about an order 
of magnitude (Tsamis et al. 2002). Other PNe central stars have presented 
important flux changes in timescales of a few years (e.g., Kostiakova 2002 and 
references therein). It is possible that this mechanism might account for the 
very low values of T(Bac), ~3000 K, found by Luo & Liu (2002) in the outer 
parts of NGC 7009. 

It would be possible also to isolate regions from the stellar ionizing photons 
if the amount of material in the line of sight increases with time, even if the 
stellar ionizing radiation remains constant. It has been suggested that dust 
clouds moving tangentially to the line of sight are responsible for the 1981-1985 
and 1996-1997 fadings of V651 Mon, the central star of NGC 2346 (Costero et 
al. 1986; Kato, Nogami, &; Baba 2001). In addition to dust clouds a similar 
effect would be caused by tangentially moving gaseous clouds of higher density 
than the ambient density. 



4. Discussion 

Temperature variations higher than those predicted by chemically homogeneous 
photoionization models of constant density are definitively present in a good 
fraction of PNe. 

Chemically homogeneous photoionization models of constant density pro- 
vide a good first approximation to the temperature structure and to the chemical 
abundances only for about half of the well observed nebulae. However, for the 
other half, it is paramount to have observations of high accuracy to determine 
their temperature structure which is needed to obtain accurate abundances of 
heavy elements relative to hydrogen. 

Spatially resolved values of Tq and arc needed to decide among the various 
possibilities mentioned above. For example, if Tq decreases in the outer regions, 
the temperature variations might be due to a decrease of the ionizing flux; in 
the opposite case the temperature variations could be due to shock heating. 

A very important observational clue to stTidy this problem is provided by 
Garnctt & Dinerstein (2001) who find a large range in the difference between the 
O"*"*" abundances derived from O II and those derived from [O III], spanning from 
no difference up to a factor of six. The size of this discrepancy is anti-correlated 
with nebular surface brightness; compact, high-surface-brightness nebulae have 
the smallest discrepancies. 

It is important to determine the cause of the temperature variations to be 
able to constrain the models of stellar evolution and of galactic chemical evolu- 
tion. If the temperature fluctuations due to chemical inhomogeneitics dominate, 
the proper abundances to use are those derived from forbidden lines under the 
assumption that ~ 0.01; alternatively if the temperature fluctuations due to 
the combination of the other six mechanisms mentioned above dominate, the 
proper abundances to use are those provided by the recombination lines. 
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The C/H values derived from recombination lines in PNe are in better 
agreement with chemical evolution models of the Galaxy than the values derived 
from collisionally excited lines under the assumption that = 0.000 (Carigi 
2002). Alternatively the abundances derived from the infrared lines seem to be 
in favor of chemical inhomogeneities, but these abundances could be higher if 
large density fluctuations are present. 

The presence of abundance gradients in the Galaxy and the values of AY/AZ 
and AY/ AO derived from galactic and extragalactic H II regions indicate that 
there are significant temperature variations in H II regions and that they are not 
due to chemical inhomogeneities (see Peimbert & Peimbert 2002 for a review). 
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